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Loop quantum cosmology tries to capture the main ideas of loop quantum gravity and to apply 
them to the Universe as a whole. Two main approaches within this framework have been considered 
to date for the study of cosmological perturbations: the dressed metric approach and the deformed 
algebra approach. They both have advantages and drawbacks. In this article, we accurately compare 
their predictions. In particular, we compute the associated primordial tensor power spectra. We 
show - numerically and analytically - that the large scale behavior is similar for both approaches 
and compatible with the usual prediction of general relativity. The small scale behavior is, the 
other way round, drastically different. Most importantly, we show that in a range of wavenumbers 
explicitly calculated, both approaches do agree on predictions that, in addition, differ from standard 
general relativity and do not depend on unknown parameters. These features of the power spectrum 
at intermediate scales might constitute a universal loop quantum cosmology prediction that can 
hopefully lead to observational tests and constraints. We also present a complete analytical study 
of the background evolution for the bouncing universe that can be used for other purposes. 
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I. INTRODUCTION 

Loop Quantum Gravity (LQG) is a consistent theory 
of quantum pseudo-Riemannian geometry that builds on 
both Einstein gravity and quantum physics, without re¬ 
quiring any fundamentally new principle (like, e.g., extra¬ 
dimensions or supersymmetry). Several introductory re¬ 
views can be found in [1]. Loop Quantum Gosmology 
(LQG) is a symmetry reduced version of LQG (see [5] 
for introductions) which accounts for the basic cosmo¬ 
logical symmetries. At this stage, a fully rigorous deriva¬ 
tion of LQG from the mother theory is not yet available. 
In fact, LQG imports the main techniques of LQG in 
the cosmological sector and uses a “LQG-like” quanti¬ 
zation procedure. This so-called polymeric quantization 
relies on a kinematical Hilbert space that is different from 
the Wheeler-DeWitt one, and therefore evades the Von 
Neumann uniqueness theorem. Nonetheless, it has been 
shown to be well defined when the diffeomorphism in¬ 
variance is rigorously imposed [3]. Since there is no op¬ 
erator associated with the Ashtekar connection but only 
with its holonomy, the basic variables of LQG are the 
holonomy of the Ashtekar connection and the flux of the 
densitized triad, its conjugate momentum. The main re¬ 
sult of LQG is that the Big Bang singularity is removed 
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and replaced by a Big Bounce smooth evolution, so that 
the total energy density cannot be greater than a criti¬ 
cal energy density. Intuitively, for sharply peaked states 
of the background geometry, the Universe undergoes a 
quantum tunneling from a classical contracting solution 
to a classical expanding solution. 

At the effective level, LQG can be modeled by two 
kinds of corrections. The inverse-volume corrections [4] 
(or inverse-triad, if one relaxes the isotropy hypothesis) 
are natural cut-off functions of divergences for factors 
containing inverse powers of densitized triads, arising be¬ 
cause of spatial discreetness. The holonomy correction 
is instead associated with higher powers of the intrin¬ 
sic and extrinsic spatial curvature components, stemming 
from the appearance of holonomies of the Ashtekar con¬ 
nection. As the status of inverse-volume correction is less 
clear -in particular because of a fiducial-cell dependance- 
we only consider in this article the holonomy corrections. 

Even when dealing with holonomy corrections only, 
there are two main ways of considering the effective the¬ 
ory, leading to a lively debate within the LQG commu¬ 
nity. This study aims at comparing the predictions for 
cosmological perturbations of both approaches, setting 
the initial conditions in the same way (that is at the same 
time and with the same vacuum), which as not been done 
to date. 

The first approach has been developed in [6H8] and is 
referred to as the dressed metric approach. It relies on 
a minisuperspace strategy where the homogeneous and 
isotropic degrees of freedom as well as the inhomogo- 
neous ones (considered as perturbations) are both quan¬ 
tized. The former quantization follows the loop approach 
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whereas the latter is obtained from a Fock-like proce¬ 
dure on a quantum background. The physical inhomo¬ 
geneous degrees of freedom are given by the Mukhanov- 
Sasaki variables derived from the linearized classical con¬ 
straints. The second order Hamiltonian is promoted to 
be an operator and the quantization is performed using 
techniques suitable for the quantization of a test field 
evolving on a quantum background [S]. The Hilbert 
space is just the tensor product of a Hilbert space for 
the background degrees of freedom, with another one for 
the perturbed degrees of freedom. In the interaction pic¬ 
ture, the Schrodinger equation for the perturbations was 
demonstrated to be formally identical to the Schrodinger 
equation for the quantized perturbations evolving on a 
classical background but using a dressed metric that en¬ 
codes the quantum nature of this background. 

The second approach, that we refer to as the deformed 
algebra, focuses on the well known problem of the consis¬ 
tency of the effective theory. This basically means that 
the evolution produced by the model should be consistent 
with the theory itself. This translates into the require¬ 
ment that the Poisson bracket between two corrected 
constraints should be proportional to another constraint. 
The coefficient of proportionality being a function of the 
fundamental variables, which makes the situation slightly 
more subtle than in usual field theories dealing with sim¬ 
ple structure constants. The key point is that the clo¬ 
sure of the algebra should also be considered off-shell 
m- Interestingly, this closure consistency condition is, 
after the holonomy correction implementation, basically 
enough to determine the structure of the quantum Pois¬ 
son bracket algebra [IIHS]. An essential result is that 
the spacetime structure eventually becomes Euclidean in¬ 
stead of Lorentzian around the bounce, when the total 
energy density is larger than half the critical energy den¬ 
sity. This had been overlooked until spherically symmet¬ 
ric inhomogeneity and cosmological perturbations were 
studied in an anomaly-free way. Without inhomogene¬ 
ity, one cannot determine the signature because (i) it is 
impossible to see the relative sign between temporal and 
spatial derivatives and (ii) the relevant Poisson bracket 
trivially equals zero in homogeneous models. The signa¬ 
ture change is not a consequence of inhomogeneities, the 
latter rather being used as a test field. There are hints 
that in the present context, such an effect could really 
be interpreted as a deep signature change of space-time 
rather than a mere tachyonic instability |14j . 

In this specific study, we do not focus on a specific 
approach. Both have their advantages and drawbacks. 
The dressed metric approach certainly captures more 
quantum effects, as it deals with the full wave functions. 
But it faces a problem. In general relativity (GR), 
there is in principle an infinite number of dynamical 
laws, all written with respect to different choices of time 
coordinates. They are all equivalent one to another 
because of the symmetries of the classical theory and 
it is legitimate to pick up an arbitrary choice. In the 
dressed metric approach, one is implicitly making use 


of several such choices, referred to as a background 
gauge. The mode dynamic is then written in terms of 
coordinate-invariant combinations of metric and matter 
perturbations. Only after these steps, one obtains 
a specific dynamic for the background variables and 
perturbations, which is written in a Hamiltonian way. 
Classically, the resulting dynamic does not depend on 
the coordinate choice and the procedure is valid. But as 
some degrees of freedom are quantized here, the equa¬ 
tions are modified by quantum corrections of different 
kinds, and nothing still guarantees that the results do 
not depend on the arbitrary choices made before (that 
is, the theory may not be covariant or anomaly-free). 
What is important is the fact that the classical theory 
enjoys a strong symmetry which is often used in order 
to simplify the analysis. When one quantizes or modifies 
the theory, this symmetry must not be violated, or else 
one may obtain meaningless (gauge-dependent) results. 
When the dynamics (including dynamical equations and 
symmetries) is formulated as a constrained system, one 
gains access to powerful canonical methods by which 
the consistency of the theory can be easily analyzed. 
It is of course possible to use another formalism, but 
not to ignore the problem of potential violations of 
crucial symmetries [14] . The deformed algebra approach 
does not suffer from this problem and is certainly more 
obviously consistent. But it does suffer from other 
difficulties, namely the shape of the modifications is not 
strictly speaking entirely determined by the anomaly- 
free condition, there is a kind of tension with the 
Hojman, Kuchar and Teitelboim theorem m making 
the geometrical interpretation difficult, and the fact that 
the fields are normalized after the effective quantum 
corrections were applied to the background, leading to a 
kind of possibly artificial “re-quantization” of the theory. 


The first part of this article is devoted to analytical 
investigations of the background evolution that were al¬ 
ready known but not expressed in such a systematic way. 
This material will also be very useful for the rest of the 
study, as the shape of the primordial tensor power spec¬ 
trum depends mainly on the cosmic history. The second 
part is devoted to the calculation of the infrared and 
ultraviolet limits of the primordial tensor power spec¬ 
trum for sharply peaked states in the dressed metric ap¬ 
proach. The third part deals with the same issues in the 
deformed algebra model. In both cases, the initial con¬ 
ditions are set in the same way, in the contracting phase, 
in order to make a meaningful comparison. The fourth 
part shows the results of the numerical computations of 
the full power spectra and some universal features are 
underlined. In the conclusion, we outline the main dif¬ 
ferences and similarities between both approaches before 
giving some perspectives towards observational tests and 
constraints. 
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II. BACKGROUND EVOLUTION: 

ANALYTICAL SOLUTIONS 

In this section we study the background evolution at 
the effective level. Although this has already been stud¬ 
ied (see i), our purpose here is to provide analytic so¬ 
lutions which are accurate approximations of the cosmic 
history over different regions. Our scope is twofold. First, 
this can give further insights on the effective dynamics of 
the background, potentially useful for further investiga¬ 
tions. Second, these analytic results are developed in the 
scope of the forthcoming investigation of tensor pertur¬ 
bations since their equation of motion obviously involves 
background quantities such as the scale factor and the to¬ 
tal energy density. Here we focus on the most probable 
dynamics as in Haim- 

A. Overview of the background dynamics 


There are two time scales involved in this system of 
equations. One is given by 1/m and corresponds to the 
classical evolution of the field. The other time scale is 
1 / y/Gpc and corresponds to the quantum regime of the 
evolution. Modulo a numerical factor, relevant for the 
following calculations, the ratio of these two time-scales 
is 

TYl 

If we assume T ^ 1, and start with a negative Hubble pa¬ 
rameter (contracting universe), the background dynam¬ 
ics splits into three subsequent phases: 

(i) Pre-bounce contracting phase, 

(ii) Bouncing phase, 

(hi) Slow-roll inflation. 


The background evolution of the quantum universe 
is described using the effective, semiclassical dynam¬ 
ics, as derived in loop quantum cosmology with holon- 
omy corrections. In this article, the background ge¬ 
ometry is described by the unperturbed metric tensor 
g = — dt (g) dt -I- a^(5ydx‘ ® dxb where a is the scale fac¬ 
tor. Dots denote derivatives with respect to the cosmic 
time, a = ^, and primes denote derivatives with re¬ 
spect to conformal time, related to the cosmic time by 
dt = adt]. The content of the universe is modelled by a 
single massive scalar field, </, with a quadratic potential, 
V {(/)) = In order to characterize the field evo¬ 

lution we use two dynamical parameters, the potential 
energy parameter, x, and the kinetic energy parameter, 
y, defined by 


m(j) 


y = 




( 1 ) 


where is the critical density, i.e. the maximum value 
of the total energy density that can be express as p = 
Pc ■ The modified Friedmann equation, as pre¬ 

dicted in LQC from the Hamiltonian constraint and the 
Hamilton equations, is 


= 



( 2 ) 


In each phase, it is possible to get analytical expres¬ 
sions for all the background variables. Note that the 
value of the inflaton mass preferred by Cosmic Microwave 
Background (CMB) observations is m ~ 10“®mp|. Fur¬ 
thermore, calculations of the black hole entropy suggests 
Pc = 0.4Imp,, leading to F ~ 2 x 10“^. Therefore, assert¬ 
ing F <C 1 is not a strong assumption at all. 


B. Initial conditions 


The initial conditions {oo, a;o: l/o} are set in the remote 
past, when Hg < 0 and 



( 6 ) 


The subscript ‘0’ means that the variables are evaluated 
at t = 0. The condition ([^ ensures that initially the 
dynamic is not dominated by the amplification due to 
the term ‘3iF’ in (|^. We often use polar coordinates for 
X and y. 


^ (^) = y^ 'sin {mt + do), 
y (^) = "v/^cos {mt + do) • 


( 7 ) 


where H = ^ is the Hubble parameter. The Klein- 
Gordon equation for the scalar field is 


The initial value of the energy density is specified with 
the two numbers a and do: 


(j) + 3H(p + m^4> = 0 . 


Equations ([^ and ^ are recast into 

( H = -SirCp^y^ (l - 2x^ - 2y^) , 
lx = my, 

I y = —3Hy — mx. 


(3) Ff sin(2do) r^^ 

V a \ 4a j 

For a given a ^ 1, such that § is valid, there is a one- 
to-one correspondence between the family of solutions to 
('^) Q and the interval {do| 0 < do < 27r}. The choice for this 
parametrization is clarified in the next section. 
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C. The pre-bounce classical contracting phase 


As long as ^ holds for p{t), the system Q can be 
solved analytically. In the third line of Q, the term 
^iHy' can be neglected, compared to mx\ as their ratio 
is of order 0{1/a) initially. Then, x and y behave simply 
as the phase variables of the harmonic oscillator {i.e. Q 
with constant amplitude). The solution for y can be in¬ 
jected into the equation for H in @ where one neglects 
‘—2x^ — 2y^’ in comparison to unity in the bracket. The 
Hubble parameter is replaced by its expression in terms 
of the energy density where the correction 1 

is neglected. After these replacements, one is left with 
a first order differential equation over p(t) which can be 
integrated into 

^ ^ [mt-b 5sin(2TOt-b26»o)] I . (9) 


This solution exhibits an oscillatory behavior due to the 
sine function in the denominator. The oscillations have 
a period of order 1 /m, much smaller than the time scale 
of the growth, ajm. Moreover, their amplitude is also 
smaller than the averaged amplitude \/piPc by a factor 
a. When these small and fast oscillations are neglected, 
the Hubble parameter can be expressed as 


H(t) = Ho (1 + l^ot)' 


( 10 ) 


where the initial Hubble parameter is Hg = —m/(3a). 
With the parametrization Q, solutions with the same a 
but different 0 o’s are all corresponding to the same aver¬ 
aged behavior (there is only a phase difference between 
them). From (10), the scale factor can be computed as a 
function of cosmic time, and as a function of conformal 
time after another integration. As the value of the ini¬ 
tial conformal time r/o be set arbitrarily, we choose 
Po = 2/(i7oao)- With such a choice, the expression for 
the scale factor simply reads 


a (p) = XoP^ with Aq = 


alHl 


( 11 ) 


so that the expression of the comoving Hubble radius 
during the contracting phase is 

aH{p) = (12) 

V 

This is the same behavior as with a universe filled with 
dust-like matter. When H ~ —mjZ^ the amplification 
term ‘3iF’ in (|^ becomes dominant. It corresponds to 
the end of the pre-bounce contracting phase and the start 
of the bouncing phase. The contracting phase ends when 
Pa = r^p,,, so at this stage there is no significant quantum 
effects. 


D. The bouncing phase 


small and fast oscillations of the field are neglected, the 
fractions of potential and kinetic energy are given by 


a:A = rsin0A and yA = rcos0A, (13) 


with 0A = 2 {a — 1) + 9o. The Hubble parameter keeps 
increasing (in modulus) until it reaches a maximum, 
= -y247rGp7/6. The inverse of has the di¬ 
mension of a time and gives an estimate of the time scale 
of this amplification. As a first analysis, in the second 
equation of the system ([^, the time derivative can be 
replaced by a factor Then, we find that the ratio 

between the fraction of potential and kinetic energy is of 
order ~ bT, and therefore very small in comparison to 
unity. This suggests that at the start of the bouncing 
phase, the kinetic energy parameter grows very quickly, 
while the fraction of potential energy remains of order 
~ r. When the kinetic energy is dominant, the system 
of equations reduces to 

iy = V247rGp„ j/Vl “ 2 /^) (-^ 4 ^ 

\x = my, 


which can be solved analytically. The solutions to (14) 
shall be valid as long as the kinetic energy dominates 
over the potential energy. In particular, they are valid 
at the bounce when the energy density reaches p^, or 
equivalently when y {t^) = 1. For the time t-g, at which 
the bounce occurs, one hnds ts = + m|cos 6 >A| ’ 

The fractions of kinetic and potential energy during 
the bouncing phase can be expressed as 


y (t) = [ 1 -b 247rC?p„(t - te)^] , (15a) 

X (t) = Xg + eFarcsinh ^-\/247rGpc (t — tg)^ , (15b) 


where e = sgn(cos 0 A), and the value of the potential 
energy parameter at the bounce is given by 


a;B = a^A — ehln (^F |cos0a|) • (16) 


The case cos 0 a "C 1 may appear problematic. Actually 
it corresponds to a different evolution of the background, 
with a phase of deflation before the bounce. Here we fo¬ 
cus on cases -statistically much more frequent and there¬ 
fore relevant for phenomenology m- where a sufficiently 
long phase of inflation is achieved. During the bouncing 
phase, the Hubble parameter and the scale factor take 
on a very simple form. The scale factor is related to the 
kinetic energy parameter by a = aelyl^s. Consequently, 
the expression for the scale factor at tA is 


= Og |r COS0A 


(17) 


Using (12), we can find the conformal time Pa that cor¬ 
responds to Then, we can use 0 and 0 in order 
to write Pa in terms of Aq. We get 


Let us define tA, the time such that H (^a) = —mf3. 
One finds ^a = 2 (a — 1) /m. Moreover, at tA, if the 



(18) 
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After the bounce, the fraction of potential energy in¬ 
creases. Meanwhile, the fraction of kinetic energy de¬ 
creases and eventually becomes smaller than the fraction 
of potential energy. This corresponds to the start of slow- 
roll inflation. 


E. The classical slow-roll inflation 


The total energy den sity p =_Pa + y^), with x (t) 
and y (t) given by (15a) and , reaches a minimum 
at time t,. According to these analytical expressions the 
total energy density increases for t > Obviously, 
this is irrelevant in an expanding universe without en¬ 
ergy sources: the total energy density must always de¬ 
crease. The time t, can be computed analytically by 
solving p(ti) = 0. One gets + (f/m), where / 

is expressed in terms of the Lambert W function (de¬ 
fined as the solution to z = and Xg is given 

by 


f = 



with 



(19) 


In general, / is of order 0{1). For instance, when 
cos^A = 1 and T = 2 x 10“^, one gets / ~ 0.18. At t,, 
the fr action of potential energy is calculated with (15a), 
(|15b|) and (fl^. We find 


X, = j;A-2erin(^iry''b^^^p5"y (20) 

Shortly after t, (in a time of order l/(r7T,lnr)), one can 
show that the fraction of kinetic energy ends up being 
almost constant. One then has y, = —eT and the slow- 
roll conditions are fulfilled. Actually, for the quadratic 
potential it is enough to check that = —H/H^ is small 
in comparison to unity for the slow-roll conditions to be 
valid. We find 


Eh 



( 21 ) 


which is generally a small number. For cos 9j^ = 1 and 
F = 2 X 10“^ one gets Ch — 0.003. Slow-roll inflation can 
start, the system of equation reduces to 


y = -^r, 

X = my, 


( 22 ) 


and the Hubble parameter becomes 

H{t) = Hi\l - e^m{t - ti)|. 


(23) 


where H, = y/STTGpJ^\Xi\. We can also use (23) to 
compute the scale factor at t, with |?/i| = F. We get 


a, = ObF 3. 


(24) 


Note that at the start of slow-roll inflation the total en¬ 
ergy density is smaller than the critical energy density 
by a factor F^. Therefore, when slow-roll inflation starts, 
the universe is already classical (since quantum correc¬ 
tions are negligible). 

We stress that all the analytical approximations de¬ 
rived above have been checked against numerical integra¬ 
tions of equation (|^. This has been done for each one of 
the three subsequent phases as well as for the matching 
between them. 


III. POWER SPECTRUM IN THE DRESSED 
METRIC APPROACH 


A. Preliminaries on the dressed metric approach 


The dressed metric approach for both scalar and tensor 
cosmological perturbations in LQC has been developed 
in [MS]. Focusing on the tensor modes, the primordial 
power spectrum at the end of inflation is defined in terms 
of the mode functions of the Mukhanov-Sasaki variables, 
denoted Vk, a^ 


V^k) 


32Gfc^ Vkiv^) 

TT a(77e) 


(25) 


with 77, standing for the end of inflation. 

It is worth mentioning that the precise knowledge of 
7^3 is not mandatory for the derivation of the primordial 
power spectrum in both the infrared (IR) and ultraviolet 
(UV) limits. For the IR limit, this is because infrared 
modes are (by definition) mainly amplified during 
the contraction and the contribution of inflation is 
suppressed as compared to the previous phases. In the 
UV, this is because we focus on modes that crossed the 
horizon during inflation, so that their amplitude has 
remained constant after a few e-folds. 


In order to obtain the power spectrum, one has to solve 
the equation of motion for the mode functions, Vk{r]), 
with given initial conditions. In conformal time, this 
equation takes the form of a Schrodinger equation 

Xkiv) + (^k^ - Mv) = 0, (26) 

where d is a dressed scale factor and (.) refers to the quan¬ 
tum expectation value on background states. This takes 
into account the width of the background wave function 
and has a priori no reason to be equal to the scale factor, 
a{t), solution to the modified Friedmann equation (cor¬ 
responding to the scale factor traced by the peak of the 


^ This model is parity invariant and the two helicity states of the 
tensor mode are equally amplified. The summation over the 
helicity states is implicitly done in our definition of the primordial 
power spectrum. 






















6 


sharply peaked wave function). However, it is argued in 
[ 8 ] that for sharply peaked background states, the dressed 
effective potential term, (a") / (a), is very well approxi¬ 
mated by its peaked value, a"/a, from the bounce up to 
the entire expanding phase. We expect this approxima¬ 
tion to be valid from the bounce down to the classical 
contracting phase since this also corresponds to a more 
and more classical universe when going backward in time 
from the bounce. With this approximation, (26 1 becomes 


v'av)+[k^--]vk{v) = 0, 


(27) 


where the scale factor is now solution to the modihed 
Friedmann equation, and the analytical results derived 
in Sec. |ll]can be used for the background variables. 


B. Calculation of the IR limit 


1. Definition of the IR regime 

The IR limit of the primordial power spectrum is ob¬ 
tained by considering the modes which stopped oscillat¬ 
ing with time and were frozen during the pre-bounce con¬ 
tracting phase. The freezing of a mode happens when its 
wavenumber becomes smaller than the effective poten¬ 
tial, a" ja. With the analytical expressions given in 
the previous section, one finds that during the contrac¬ 
tion. 


a 

a 


1 


2 ■ 


(28) 


Thus, an infrared mode with a wavenumber k crosses 
the effective potential at a conformal time | 77 fc| = jk. 
Its amplitude is frozen from that time up to the end 
of inflation, as k'^ remains smaller than a"fa. Since 
—00 < T] < r]A (with r]A < 0), the modes that crossed the 
potential during the contracting phase are in the range 
0 < fc < fciR, with ki^ defined by the mode that crossed 
the effective potential at the beginning of the bouncing 
phase. With (HU), ([g, ^ and ( [T^ we find 


j Ob f m?^/24frGpf 

^ 1 |cos 0 ,| 


1/3 


(29) 


The IR limit stands for the modes such that k ki, 


where ak and Pk are two constants to be determined. 
The value of ry* can be conveniently set by requiring the 
term proportional to ak to be solely decaying, and the 
term proportional to rjk to be solely growing. During the 
contracting phase, the term proportional to ak is clearly 
decaying since 0 ( 77 ) is decreasing. A convenient choice 
for is such that the term proportional to fik niust be 
solely growing. Since 0 ( 77 ) = XoT]^, the term proportional 
to Pk has a time dependance ~ 77 ^( 77 “^ — in which 

the part proportional to 77 ^/ 77 ^ is decaying (77 < 0 ) and 
we send 77 * to (— 00 ) to remove ii[^ 

From (30) and (25), the expression of the IR limit of 


the spectrum reads 


'PT(fc)™ = \ak + Pkli'n^)f : 


(31) 


where /(t/^,) is the integral defined as 


(32) 


The calculation of the IR limit proceeds in two steps. 
First, we compute ak and Pk by matching (30) to a set 


of solutions defined in the contracting phase. As we 
shall see, this determines the scale dependence of the 
primordial power spectrum in the infrared regime. The 
second step is the calculation using the analytical 

solutions for the background, obtained in Sec. |ll] This 
second step sets the amplitude of the power spectrum. 
The expression of the primordial spectrum is finally 
obtained by gathering the expressions of ak, Pk and 
HVe)- 


In order to derive the expressions of ak and Pk, the 
approximate solution given in (301 (valid in the IR only 
but from 77 ^ to 77 ,) has to be matched with a set of solu¬ 
tions to the equation (27), during the contracting phase. 
With a” fa = ‘Ifrp, this set of solutions corresponds to 
the linear combinations of the Hankel functions of order 
V = 3/2: 


= V-kv [Afci73/2(-/c?7) + BkHl.^{-kT])] , (33) 

where the superscript ‘C’ recalls that ( [M| ) is valid only 
during the contracting phase. In order to specify Ak 
and Bk we match ( |33| ) with the Minkowski vacuum in 
the remote past, i.e. Vk{r] —>• — 00 ) = This 

requirement leads to 


Ak 



Bk = 0 , 


(34) 


2. Primordial power spectrum in the IR regime 


For infrared modes, from the potential crossing rjk to 
the end of inflation 77 ^, the solution to the equation of 
motion (27) is therefore well approximated by 


=akair]) + Pka{r]) f -^^^+0{(k/k,^f), (30) 
Jm “ ) 


^ During the contracting phase, the identification of the growing 
and decaying modes differs from that identification during in¬ 
flation. Because the Universe is expanding during inflation, the 
term akCL{r}) is solely growing (while it is solely decaying dur¬ 
ing contraction). Then the term dr ]'can be made 

solely decaying in an inflationary universe by setting rji, = r]e 
(while it is made solely growing during contraction by setting 
?7* ^ -00). 
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up to a phase which is has no importance her^ A set 
of solutions valid in the range — oo < 77 < 77 ^, and corre¬ 
sponding to the Minkowski vacuum, is thus 




(35) 


Since 77 ^ <C 77 ^ for infrared modes, the IR limit of (35) 


has to coincide with (30) in the interval Tjk r] < rjj^. 


At a given rj in this interval, we calculate the asymptotic 
limit of the Hankel function when fc —)■ 0. This leads to 

= (=>«> 

The term of order has a time dependence given 

by 0 ( 77 ) oc 77 ^, and corresponds to the term proportional 
to ak in (30). 


Eventually, we have to match (36) with the explicit 
expression of (30) that one obtains with 0 ( 77 ) = and 
77 * = - 00 : 


vTiv) = Oik\rf - 


Pk 

3Ao?7 


(37) 


By comparing (|37[) with (36), one finds 


ak=0{k^/'^) and /3fe = (3z/V?)AoA:”^/^. (38) 

For infrared modes the contribution of Uk is negligible, 
so that (31) simplifies to 


^T(fc)“ = ^A2|/(77jr 


(39) 


Therefore, in the IR limit we expect the power spectrum 
to be scale invariant (at least at the order of validity of 
our approximations). 

The amplitude of the power spectrum in the IR regime 
is obtained by evaluating the integral I{r]^) = drj/a^. 
In order to do this, we split the integral into three parts 


Hv.) = -f/(77A,??i) -b/(77i,77e). 


(40) 


y = x/m, the integrand is proportional to x and the inte¬ 
gral itself is proportional to the difference \Xi — Xf^\ (which 
is given in (20)). Eventually, one gets 


= - 


m 


1 


I8An |cos6*a 


In . (42) 


The last part of the integral corresponds to the slow- 
roll inflation as obtained from a massive scalar field. The 
calculations are well known in this case, leading to 




1 


I 


3a3i7, 3a377, 


[l + 0(e„)], (43) 


where Ch = —H/H^ is the slow-roll parameter which re¬ 
mains small in comparison to unity (except in the neigh¬ 
bourhood of te). It will be neglected in the forthcom¬ 
ing calculations. During slow-roll inflation, the Hub¬ 
ble parameter decreases linearly with cosmic time while 
the scale factor grows exponentially. The second term, 
l/(afi7e), can be safely neglected as it is suppressed by 
a factor ^ exp(—3iV,), where denotes the number of 
e-folds from 77 , to 77 ,. This also means that the detailed 
dynamic of inflation is not needed here, since its contri¬ 
bution is rapidly negligible after a few e-folds. With the 
expressions of a, and 77, given in ( [^ , I{rii,r]^) evaluates 
to 


I = 


r 


12 v^Ao |a:iCos6>A 


(44) 


with Xi given in ( 20 ). 


Gathering the results (41), (42) and (44), the integral 
7 ( 77 ,) can be written as 7 ( 77 ^ = (1 -|-1 -|- J'), so that 


the IR limit of the power spectrum (39) reads 


r^{ky^ = ^m^\i+i+j\^, 

yTT 


(45) 


where 


The first part corresponds to the contracting phase, the 
second part corresponds to the bouncing phase, and the 
last part gives the contribution of the inflationary phase. 
With 0 ( 77 ) = Xgrj^ during the contracting phase, and re¬ 
calling that 77 a = —[ 6 /( 77 iAo)]^/^, the first part of the 
integral is easy to compute: 

/(-oo,77a) = ^. (41) 

The second part of the integral is first written in cos¬ 
mic time, 7 ( 77 A, 77 i) = During the bounc- 

_ 1 

ing phase we have found that a = ObIt/I 3. Then, with 


3 This also fits with the appropriate Wronskian condition as re¬ 
quired for the quantization a la Fock of the tensor perturbations 
field. 



(46) 

icos0a|^"1^^v / 


(47) 

2 X; COS 0A 


In general, J' is much smaller than I, suggesting that the 
contribution to the amplitude of the spectrum in the IR 
that corresponds to inflation is negligible (for instance 
with cos 0 A = 1 and T = 2 x 10 “^, one gets JjX ^ 0 . 002 ). 

The scale-invariance of the IR limit of the spectrum is 
a direct consequence of the fact that the infrared modes 
crossed the effective potential, a"/a, during the contract¬ 
ing phase whose dynamics is equivalent to that of a dust¬ 
like matter dominated era. No further assumption on the 
detailed dynamics of the bounce is needed to get the scale 
invariance (though the detailed dynamics is needed to get 
the amplitude of the power spectrum). The amplitude 































only depends on three parameters: the critical energy 
density, the mass of the scalar field, and the phase 0 a, 
between x = m(j)/y/2p^' and y = (j)/ y/2p^ at the start of 
the bouncing phase. The first two parameters are funda¬ 
mental. The phase 0 a depends on Og which is a contingent 
parameter whose value sets the initial conditions (see m 
for a more detailed discussion). The case cos0 a 1 
may appear problematic (as it would lead to a divergent 
power spectrum), however in this case the dynamic of the 
background would be different (with deflation before the 
bounce) and our analytical results would not be valid. 


C. Calculation of the UV limit 

1. Definition of the UV regime 

By definition, the ultraviolet modes have remained 
well inside the Hubble radius until the phase of slow-roll 
inflation. They are insensitive to the background cur¬ 
vature during the contracting and the bouncing phase. 
The effective potential a"/a can be written in terms of 
the Hubble parameter and its time derivative as a"/a = 

tff (^H + 2H^^. During the bounce, this expression be- 
cornes 


SirGp^ 


aly'^ (V - 1 ). 


(48) 


It is clear in (481 that the effective potential reaches its 


maximum at the bo unce, w hen y = 1. This feature sets 
a scale, = max a/ a" ja \ which evaluates to 


fcuv = Oba/SttG^. 


(49) 


commonly used in slow-roll inflation. The power spec¬ 
trum in the UV regime is then given by the standard 
red-tilted power spectrum of slow-roll inflation (see HS¬ 
US]), 


dlnfc 


16G 

TT 


[1 


— —Sch, 


2eH(2C + l)], (51a) 
(51b) 


where H is the Hubble parameter evaluated when k = 
aH, and C ~ —0.73. At the order of validity of our ap¬ 
proximation (r ^ 1) and neglecting Ch in the amplitude, 
these expressions become 


V^ikf'' 

dln/c 


16G 2 
=- m 

TT 


Xi 

r 



2 


(52a) 

(52b) 


where T = mj\/2AfdTf)f and x, is given by (20). This 
prediction of a slightly red tilted spectrum matches the 
standard inflationary model. The amplitude scales with 
rrif and depends on the critical energy density in a non¬ 
trivial way. With cos 0 a = 1 for simplicity, one gets 
'P’rikY'^ cx Ilf {mjg/Gpf ). Moreover, with the stan¬ 
dard value r = 2 X 10“^ we find that the spectral index at 
the start of inflation, given by (52b), is Ut — 1 — —0.007. 


IV. POWER SPECTRUM IN THE DEFORMED 
ALGEBRA APPROACH 


A. The deformed algebra approach 


All modes with a wavenumber larger than fcuv crossed 
the potential during slow-roll inflation. The UV limit 
of the power spectrum is defined by the modes with a 
wavenumber k ^ fcuv 

2. Primordial power spectrum in the UV regime 


The calculations presented above can be extended to 
the case of the tensor power spectrum in the deformed al¬ 
gebra approach [TMH]. The equation of motion for the 
mode functions of the Mukhanov-Sasaki variables also 
takes the form of a Schrodinger equation. However, the 
frequency term is time-dependent and the effective po¬ 
tential is different: 


In the dressed metric approach, the calculation of the 
UV limit of the power spectrum is straightforward. As 
during the bouncing phase the mode functions do not feel 
the curvature of space-time, they are well approximated 
by 


vriv) = for p < ry, (50) 

Once the Universe enters inflation, the term a"/a cannot 
be neglected anymore and behaves as (2 -|- 2>e^)/p^. The 
mode functions are now given by a linear combination of 
the Hankel functions of order 3/2-|-eH. At this stage, the 
derivation of the primordial spectrum is simple: we have 
to match the Minkowski vacuum (well defined within the 
Hubble radius for k 3> k^^) with the mode functions 


(^k^-^'^Vk{p)=A: (53) 

where 

fl = 1-2-^ and (54) 

Pc VO 


The region with H > 0, corresponding to p < Pc/2, is 
Lorentzian whereas the region with H < 0, corresponding 
to p > Pc/2, is Euclidean. Here, the mode functions 
are related to the amplitude of the tenso r mode s of the 
metric perturbation, via Vk = Zr^hk/ f A2 ttG , so that 
the power spectrum is now defined as 


V^{k) 


32Gk^ Vkivf 

TT 


(55) 
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where 77 , denotes the conformal time at the end of slow- 
roll inflation. Actually, this definition is equivalent to 
(25) because during slow-inflation z,, — a (as a conse¬ 
quence of Pi <C Pc). 


the integral. Indeed, recalling that before ?7 a the energy 
density remains smaller than we have 

/a(-oo,77A) =-2 / —^<-2T^I{-oo,riA)- (60) 

7_oo Pc 


B. Calculation of the IR limit 


The same holds for In(rii,r]^), with p, instead of Pa. The 
remaining part of the integral is 


The IR limit is defined exactly in the same way as in 
the dressed metric approach. From the expression of a 
and p as functions of conformal time, one easily notices 
that — z'r^jzT ~ — 2/rf + 0(T‘^/rf). In the con¬ 

tracting phase, there is therefore no noticeable difference 
between the deformed algebra and the dressed metric ap¬ 
proaches. The IR limit still corresponds to modes with 
k /ciR, where fem is given in (29). 

The calculation of the IR limit of the spectrum pro¬ 
ceeds in the same way as for the dressed metric approach. 
We first write the approximate solution to (53) in the in¬ 
frared regime. 


dn(pA,Pi) =-2 / (61) 

JriA Pc 

During the bouncing phase, p = and a = aB\y\~^ so 
when we switch to cosmic time, the integral becomes 


In{VA,vd = -i‘^/al) \y{t)\^dt. 


(62) 


' (a 


The last step is to express dt in terms of dy with (14). 
Then the integration can be performed analytically and 
leads to 


= ctkZTiv) + PkZTiv) 



drj' 

4(p') 


+ 0(fc2), (56) 


In{VA,Vi) 


m sin^ Pa p 2 

36Ao cos Pa 


(63) 


from which the general expression of the IR limit of the 
power spectrum directly follows: 


Vr{ky^ 


32GP 

TT 


oik + Pk 



2 


(57) 


With the definition of z^, the integral on the RHS is 
simply given by the sum of /(p,) -I- 7n(Po), with /(po) 
defined in (32), and 


J-00 Pc 


(58) 


As before, the two constants ak and Pk in (57) are ob¬ 


tained by matching the solution (56) with a set of solu¬ 


tions valid for any wavenumber k during the contracting 
phase. During the contracting phase the difference be¬ 
tween (D/c^ — z"/zt) and (fc^ — a”/a) can be neglected. 
Consequently, the two constants ak and Pk take the same 
value as before: ak = C>(fc^/^) and Pk = {3i/V^)Xok~^/^. 
With these expressions, the IR limit of the power spec¬ 
trum becomes 


^T(fc)“ = —A2|/(pj + /n(pjr 


(59) 


Note that all the differences between the deformed alge¬ 
bra and the dressed metric approach are encoded in the 
integral IniVe)- 

Now, we will show that /1(r]^) = so that 

the contribution of InijlP) to the IR limit of the spectrum 
can be neglected. First, we split the integral into three 
parts: /a(Pe) = 7n(-oo, Pa)- f/^( pa, Pi) + dn(pi, Pe). The 
proof is straightforward for the first and the third parts of 


Therefore, at order ©(F^), we predict no difference for 
the IR limits of the power spectra in both approaches. 
The IR limit of the power spectrum in the deformed al¬ 
gebra approach is still given by (45). 


C. Calculation of the UV limit 

The UV limit of the power spectrum in the deformed 
algebra approach has already been discussed in here 
we recall the conclusion of this previous work. Thanks to 
numerical integrations for the equation of motion as well 
as WKB based arguments, it is shown that the primordial 
power spectrum exponentially grows with the wavenum¬ 
ber fc, for large values of k. Actually, oscillations are 
still superimposed to this exponential envelope. During 
the bouncing phase, the term z"/z^ reaches a maximum 
|z"/zT|tB = 407rGpc- This means that for modes such 
that k^ ^ dOnGp^, the time-dependent frequency in the 
equation of motion, flk^ — z"/zt, is dominated by flk^ 
during most of the cosmic history prior to inflation. The 
Euclidean phase around the bounce, D < 0, leads to an 
instability in the equation of motion so the amplitude 
of tensor modes recieves a real exponential contribution, 
i.e. Vk^oD cx; exp(A: x \/l^l dp), where the integra¬ 
tion is performed over the interval Arj corresponding to 
the Euclidean phase. 


V. POWER SPECTRUM AT ALL SCALES: 
NUMERICAL RESULTS 

Deriving the primordial power spectrum for tensor 
modes at all scales requires a numerical integration of 
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both the equation of motion of the mode functions ((Izl) 
or (53) depending on the approach) and the equations 
of motion for the background (gathered in the system of 
equations @)- The numerical integration is performed 
starting in the contracting phase, when Po/Pa ^ T^. For 
the background, the initial conditions are set by choosing 
the value of 9g, the initial phase between the share of 
potential energy and kinetic energy in the total energy 
density. For the perturbations, the initial conditions are 
set during the contracting phase when the modes are well 
inside the horizon. The initial state of the perturbation 
can then be identified with the usual Minkowski vacuum. 


The detailed dynamics of the background {e.g. the 
value of X at the bounce, or the number of e-folds during 
inflation) and subsequently the detailed shape of VT{k), 
are fully determined by two types of parameters: the 
mass of the scalar field and the critical energy density on 
one hand, and the phase 9o on the other hand. 

The mass of the scalar field and the critical energy den¬ 
sity can be seen as fundamental constants of the model. 
Though their values are not known, some particular val¬ 
ues are favored by CMB observations and theoretical con¬ 
siderations. Even if the details of the calculation using 
the minimal area gap of LQG still need clarification, some 
dimensional arguments lead to believe that the value of 
Pc should not be far from the Planck scale. Note that p^ 
is the only parameter linked to LQG (via its dependence 
on the Immirzi parameter, 7 ). The value commonly ac¬ 
cepted is Pc = O.dlTOpi, and we shall use it as the stan¬ 
dard choice in our numerical simulation. The value of 
the mass of the scalar field, as deduced from the GMB 
observations, is generally chosen to be m ~ 1.2 x 10 “®mpi 

m- 

The parameter 9o has a different status since it is 
totally contingent. Its value can vary between 0 and 
27r (actually the range 0 < 6*0 < tt is enough as the 
equations remain unchanged under the transformation 
9o ^ 9o tt). As underlined in HZ], most of the values 
of 9o lead to a universe with a phase of inflation shortly 
after the bounce (and no deflation before the bounce). 
We have restricted ourselves to this kind of solutions 
since they are the most probable, and the more in line 
with our current knowledge of the cosmic history (be¬ 
lieved to have underwent a phase of primordial inflation). 

Qualitatively, we can already anticipate the global 
shape of the primordial power spectrum. Irrespective 
of any approach, its shape is driven by the background 
evolution through the functions a and 12 (in the de¬ 
formed algebra approach) and their time derivatives. 
Our analysis is restricted to the wide range of cosmic 
histories that split into three main eras: a classical 
(dust-like) contracting phase, a bouncing phase where 
quantum effects are significant, and a classical inflation¬ 
ary phase. We anticipate the shape of the primordial 
power spectrum to be qualitatively unaffected by the 
values of m, p^, and 9o. (Obviously, the precise values 


of these parameters will affect the scales and amplitudes 
involved in the spectrum at a quantitative level). We 
can also anticipate three regimes in the power spectra, 
corresponding to: the modes that have left the horizon 
during the contracting phase (large scales); the modes 
that have left the horizon during the bouncing phase 
(intermediate scales); the modes that have remained 
within the horizon until the start of the inflationary 
phase (small scales). For the large and small scales, we 
should recover the IR and UV limit derived analytically 
in the previous sections. 

In the next three sections we present the primordial 
power spectra obtained within each approach. We study 
the influence of the three parameters, m, Pc, and 9q. For 
each varying parameter, we present the primordial power 
spectra as predicted by each approach, thus facilitating 
the comparison. 

We use Planck units hereafter, with the following def¬ 
inition of the Planck mass, mpi = l/v^. For simplicity, 
we normalize the scale factor at the time of the bounce, 
setting Ob = 1. The power spectra are depicted as func¬ 
tions of the comoving wavenumber k. 


A. Varying the mass of the scalar field 

The primordial power spectrum for different values of 
the mass of the scalar field is presented in Fig. The 
upper panel corresponds to the dressed metric approach 
and the lower panel to the deformed algebra approach. 
The mass takes three different values: m = 10“^mpi (tri¬ 
angles), m = 10 “^-^mpi (open disks), and m = 10 “^mpi 
(black disks). For numerical convenience, these values 
are larger than the preferred value. However, the results 
can be extrapolated and the associated phenomenology 
shall be studied with values closer to 10 “®mp, [28]. 
The critical energy density is set equal to O.dlrripi and 
COS0A — 1- 


In the dressed metric approach (upper panel of Fig.[^, 
there are three regimes in the primordial power spectrum: 

(i) At the largest scales, for k < fcip with 

1/3 


kin — 


1 


/ m^-\/247rGpc 

3V^ V |cos6»a| 


the power spectrum is scale invariant in agreement 
with the analytical calculations of Sec. m This cor¬ 
responds to modes that were amplihed mainly dur¬ 
ing the classical contracting phase. In Fig. the 
scale corresponding to fcip is depicted with vertical 
dotted lines. At this scale, there is a transition in 
the numerical results that is in perfect agreement 
with the analytical formula (especially its de¬ 
pendence). Moreover, it is clear on the figure that 
the numerical IR limit of the spectru m b ehaves as 
m?, again in perfect agreement with (45). 
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(ii) For intermediate scales, such that < k < k^^, 
with 

fcuv = y/SirGp^ 

the amplitude of the power spectrum is oscillating. 
This part of the spectrum corresponds to modes 
that were amplified during the bouncing phase. The 
first peak corresponds to a maximum of the power 
spectrum, that reaches 100 x approximately. 
Then, the amplitude of the oscillations is damped 
for increasing values of the wavenumber. Intuitively, 
these oscillations can be understood as due to quasi¬ 
bound states in the effective Schrodinger equation. 
The second transition scale, fcuv, is depicted in 
Fig.[T] as a vertical dashed line and is in agreement 
with the transition scale found numerically. 

(iii) At the smallest scales, k > fcuv, the power spectrum 
is a power law with a slightly red spectral index, 
just as predicted by the standard inflationary 
paradigm. This part of the power spectrum cor¬ 
responds to modes that have remained inside the 
horizon until the start of the inflationary phase. 
The numerical results are in agreement with the 
UV limit derived analytically in Sec. |III[ see Fig. 

In the deformed algebra approach the primordial ten¬ 
sor power spectrum also features three different regimes 
(see the lower panel of Fig. [^. The first two regions 
(j.e. the large scales, k < k^^, and the intermediate 
scales, /ciR < k < k^^), are almost identical to the power 
spectrum derived in the dressed metric approach. The 
scale-dependence of the power spectrum and the tran¬ 
sition scales are the same. This is because the impact 
of n is subdominant for these modes. However, within 
these two regions the numerical results suggest that for 
k < fcuv the amplitude of the spectrum in the deformed 
algebra approach is slightly smaller than in the dressed 
metric approach (by less than a factor 2). (This feature 
could not be explained by our analytics.) 

At smaller scales, k > k^^, the primordial power spec¬ 
trum in the deformed algebra approach strongly differs 
from the one predicted by the dressed metric approach. 
As alrea dy s uggested by our analytical considerations 
(see Sec. IV), the power spectrum is exponentially in¬ 
creasing with the wavenumber (as a result of the insta¬ 
bility generated by Q which is negative-valued around 
the bounce), with superimposed oscillations. Note that 
the numerical results confirm once again that the scale 
defining the transition between the intermediate scales 
(oscillations) and the large scales (exponential growth) 
does not depend on to. 

The UV behavior of the spectrum clearly raises ques¬ 
tions. The first one is related to the fundamentally trans- 
Planckian nature of these modes. As demonstrated in [6], 
this is not a problem when considering the appropriate 
length operator in loop quantum gravity. A more seri¬ 
ous issue is related to the use of the perturbation theory 


Pr(k) 



Pr(k) 


•, o 

- o** 



Figure 1. Primordial power spectra for tensor modes in the 
dressed metric approach (upper panel) and in the deformed 
algebra approach (lower panel) for different values of the mass 
of the scalar field. The critical energy density is pc = O.llmpi 
and cos 0 a — 1- The mass of the scalar field takes three values; 
m = 10~®mpi (triangles), m = 10~^'®mpi (open disks), and 
m = 10“^mpi (black disks). The dashed vertical line at large 
k, corresponds to fcuv (which does not depend on m). The 
dotted vertical lines at smaller k correspond to fcm (which 
scales as 


when the spectrum increases exponentially. Obviously, 
backreaction should be taken into account in this regime 
and the results shown here are not fully reliable anymore. 
They just give a general trend and not the accurate shape 
of the spectrum. However, we believe that this is basi¬ 
cally enough for the phenomenological purposes we are 
interested in. The most interesting region, that is the 
oscillatory one, is under control and the Cg CMB spec¬ 
trum can be safely calculated [55] • If the observational 
window of wavenumber was to fall on the exponentially 
rising part, this would anyway lead to a situation incom¬ 
patible with data (as the tensor to scalar ratio is small). 
The perturbation theory breaks at a level where tensor 
modes are anyway excluded by current data. 
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B. Varying the critical energy density 

The critical energy density depends on the Immirzi 
parameter, a fundamental parameter in LQG, whose 
value is traditionally deduced from a calculation of the 
black holes entropy. Nonetheless, it has been recently ar¬ 
gued [3D] that the formula for the entropy of black holes 
can be recovered, in the framework of LQG, without 
specifying the value of the Immirzi parameter. Recently, 
a quasi-local description of a black hole m was indeed 
shown to allow one to recover at the semi-classical limit 
the expected thermodynamical behaviors of a black 
hole for all values of 7 ED, assuming the existence 
of a non trivial chemical potential conjugate to the 
number of horizon punctures. A detailed microscopic 
mechanism was also put forward in |32j and [33] where 
the area degeneracy was analytically continued from 
real 7 to complex 7 and evaluated at the complex values 
7 = ±L This motivates us to consider other values for 
the critical energy density and discuss how it can affect 
the primordial tensor power spectrum. 

In Fig. we show the primordial tensor power spectra 
for different values of p^. Here, the mass of the scalar 
field is set equal to m = 10“^mpi, and cos0 a ~ 1. The 
upper panel corresponds to the dressed metric approach 
and the lower panel to the deformed algebra approach. 
The different values of the critical energy density are 

= O.OOdiTOpi (triangles), = O.OdlrUpi (open disks), 
and Pc = 0.4Imp, (black disks) which is the theoretically 
favored value. 

The global shape of the primordial power spectrum 
is recovered for both approaches, with three different re¬ 
gions. The positions of the transition scales, fcjp and fcpv, 
clearly depend on p^ irrespectively of the approach. The 
IR transition scale, /cm, mildly decreases for smaller val¬ 
ues of Pc in agreement with the analytical calculations 
that led to fcip oc (Gp^)^/®. The UV transition scale, 
is more strongly dependent on the value of the critical en¬ 
ergy density, also in agreement with the scaling derived 
analytically, k^v oc V^- 

For the dressed metric approach, a decrease of p^ 
yields a slight decrease of the amplitude of the primor¬ 
dial power spectrum at all scales. This feature is also 
suggested by the analytical results, as both formulae 
for the UV and IR limits depend on the critical energy 
density as ^ ln^(m/-\/Gp„'). 

For the deformed algebra approach, a decrease of p^ 
leads to a slight decrease of the amplitude of the spec¬ 
trum at large and intermediate scales as in the dressed 
metric approach. At smaller scales, k > /cuvi the small¬ 
est value of Pc corresponds to the fastest divergence of 
the spectrum. Analytically, we expect this divergence 
to scale as oc exp(fc dp), where the interval Arj 

corresponds to the euclidean phase. Therefore we can 
define the rate of growth of the spectrum in the UV as 

kn = 1 / ^/Wdp- 


Pr(k) 



Pr(k) 


0.1 




0.01 



Figure 2. Primordial power spectra for the tensor modes in 
the dressed metric approach (upper panel) and in the de¬ 
formed algebra approach (lower panel) for different values 
of Pc. The mass of the scalar field is m = 10“^mpi and 
COS0A — 1. The critical energy density is Pc = O.OOdlmpi 
(triangles), Pc = O.Odlmpi (open disks) and pc = O.dlmpi 
(black disks). 


This integral can be computed, leading to 

kn ~ 0.87/247rGpc. (64) 

So, small values of the critical energy density indeed cor¬ 
respond to a quicker divergence of the power spectrum 
in the UV. 


C. Dependence on Oo 

The primordial power spectra for different choices of 0o 
are shown in Fig. |3 in the dressed metric approach (up¬ 
per panel) and in the deformed algebra approach (lower 
panel). The mass of the scalar field is to = 10“^TOpi, 
and the critical energy density is Pc = 0.4lTOp,. We chose 
five values of Og, equally spaced between {■k/2 — 1 ) and 
{■k/2 -b I), ensuring that the background goes through 
a phase of inflation after the bounce. In the numerical 
simulations we have always set a = 177r/4-|-I, so that 
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Pr(k) 


0.01 




Figure 3. Primordial power spectra for the tensor modes in 
the dressed metric approach (upper panel), and in the de¬ 
formed algebra approach (lower panel). The parameter 6o 
varries from ( 7 r /2 — 1) ~ 0.18 x tt to ( 7 r /2 -f 1) ~ 0.81 x tt. The 
exact values are: plain disks for Og = (7r/2 — 1), open disks 
for 6o = (ti/2 — 1/2), plain triangles for 9o = 7 r/ 2 , open trian¬ 
gles for 9o = (' 7 I '/2 -f 1 / 2 ), and plain stars for 9o = ( 7 r /2 -|- 1 ). 
The mass of the held is m = 10“^mpi, and the critical energy 
density is = O.llmpi. 


COS0A = 0 (with 0A = 2(a — l) + 0o) corresponds to = 0 
and a is significantly larger than one. 

The main impact of Og is in the IR regime. Both the 
infrared transition scale fcip (varying as ~ |cos(0a)|~^^^) 
and the amplitude of the IR limit of the spectrum are 
significantly depending on 0g. This is true in both ap¬ 
proaches. At intermediate and smaller scales, the power 
spectra are nearly independent of the choice of Og, again 
irrespectively of the considered approach. The numerical 
results confirm that the ultraviolet transition scale fcpv 
is independent of 0g. Moreover, in the deformed algebra 
approach the growth rate of the spectrum in the UV 
appears to be independent on 0g too, in agreement with 

p). 


In order to highlight the dependence of the IR limit as a 
function of Oq, Fig. shows this limit in both approaches 
and for different choices of 0g, with m = 10“^mpi, and 


D IR. UV 



2 


Figure 4. The infrared limit (disks) and the ultraviolet limit 
(triangles) of the primordial power spectrum as a function 
of 9g. The solid black curve corresponds to the analytical 
calculations for "Pt™, see (451. The IR limit from a numerical 
simulation is displayed with open disks for the dressed metric 
approach, and black disks in the deformed algebra approach. 
The dashed black curve stands for the analytical UV limit 
in the dressed metric approach, see ( |52a[ ). The UV limit in 
the dressed metric approach, as derived from the numerics, 
corresponds to the black triangles. The mass of the scalar 
field is m = 10 “®mpi, and the critical energy density is pc = 
0.41ml. 


Pc = 0.41mp|. The solid black curve corresponds to 
the analytical calculation for V-r™, see (451. This an¬ 
alytical curve is valid for both the dressed metric and 
the deformed algebra approaches at first order in T = 
m/^/24TrGpc. The numerical derivation of the IR limit is 
displayed as open disks for the dressed metric approach, 
and as black disks in the deformed algebra approach. We 
observe a fairly good agreement between analytical and 
numerical results. Although there are some differences in 
the amplitude of the IR limil|^ the behavior as a function 
of 9g is consistent between the analytics and the numerics. 
This shows that 7^^™ strongly depends on 9g, the former 
varying by more than one order of magnitude from its 
minimal value at 9g = 7r/2 (thus giving cos 0 a = 1)? to 
its maximal value reached when Og tends to 0 or tt. 

In the restricted case of the dressed metric approach, 
the UV limit as a function of 0o is also displayed in Fig.|^ 
The dashed black curve stands for the analytical calcu¬ 
lation presented in (52a). The UV limit obtained from 
the numerical simulation is displayed with triangles. A 
good agreement is also observed here. Nonetheless, the 


^ The discrepancy is not surprising. First of all, the analytic result 
is based on some approximations for the time dependence of a 
and n. Second, the numerical evaluation fo cannot be 

exactly obtained for k = 0 since this would require to start the 
numerical integration at 77 ^ —00 which is unfeasible. We believe 
that these features are at the origin of the disagreement between 
the numerics and the analytics. 
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remaining difference between the analytical and numer¬ 
ical results certainly comes from the approximations in¬ 
volved in the determination of Xj, on which the UV limit 
depend. 

VI. CONCLUSION AND DISCUSSION 

In this work, we have compared the dressed metric 
and deformed algebra approaches to loop quantum cos¬ 
mology. In order to compare them efficiently, we have 
set the initial conditions in the same way for both ap¬ 
proaches (in the remote past of the classical contracting 
branch). This is consistent and arguably the most ob¬ 
vious choice if the word initial is to be taken literally. 
It is however fair to mention that this is not the only 
choice one could have made. As far as the dressed met¬ 
ric approach is concerned, the authors who developed the 
strategy have preferred to set the initial conditions at the 
bounce [SHE]. Then the initial state for tensor perturba¬ 
tions is given by a 4th-order WKB vacuum defined for 
k > kuv In fact, their results seem to be very similar 
to ours (for the range of scales covered by both choices 
of initial conditions). As far as the deformed algebra ap¬ 
proach is concerned, it should be underlined that it is also 
possible to set initial conditions at the surface of signa¬ 
ture change. This has been investigated in [5S] and leads 
to a different spectrum. If these issues are left for future 
considerations and if we focus on the comparison with 
similar initial conditions, several important conclusions 
can already be drawn. 

First, it is remarkable that for both approaches the IR 
limit is the same and basically agrees with the prediction 
of standard general relativity. Therefore at the largest 
scales, the primordial tensor power spectrum cannot be 
used to probe quantum gravity (at least in this setting). 

Second, there is a strong difference between the ap¬ 
proaches in the ultraviolet regime. Whereas the dressed 
metric simply leads to the slightly red-tilted power spec¬ 
trum, as predicted in standard inflationary cosmology, 
the deformed algebra leads to an exponentially increas¬ 
ing spectrum (modulated by oscillations). 

Third, at intermediate scales, a very interesting behav¬ 
ior appears. Not only because it is substantially different 
from the predictions of the standard inflationary models 


but also because both predictions are in agreement with 
each other! This region seems to exhibit a universal 
LQC effect that has been searched for during the last 
decade. In addition, the phase of the oscillations that 
appear at these intermediate scales, does not depend 
on the unknown (and fundamentally random) phase 
parameter, Og. This opens an interesting avenue in the 
perspective of testing the predictions of effective LQC. 

In the future, this work should be extended in two di¬ 
rections. One is to consider not only the tensor modes, 
that have not yet been observed, but also the well known 
scalar modes. The relevant equations have already been 
derived for the dressed metric approach but are still to 
be investigated into more details in the deformed alge¬ 
bra approach. The reason for this difficulty is related to 
divergences (at the bounce and at the change of signa¬ 
ture) that should be regularized. The difficulty is how¬ 
ever more technical than conceptual and should be solved 
soon. 

The second path to follow is naturally to go more 
deeply into the phenomenology of this comparison 
and calculate the corresponding cosmic microwave 
background Cg spectra which are already constrained by 
observations. Two main tasks will have to be pursued. 
The first is related to the number of e-folds that the 
Universe underwent since the bounce. This number 
depends, among other parameters, on 0g and on the 
reheating temperature. Once the number of e-folds since 
the bounce will be specified, the range of wavenumbers 
considered in this study that falls within our observable 
window will be completely determined. The second 
important task is to investigate how the oscillations at 
intermediate scales shall be washed out by the transfer 
phenomena that occur between the end of inflation and 
the decoupling. 
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